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Abstract — We propose an experimental study of 
parametric resonance between P-waves and S-waves, 
which can be used to describe various nonlinear 
phenomena qualitatively and to obtain bifurcation 
diagrams quantitatively. It is based on an electronic 
circuit and is easy to design. We show that it is a good 
simulation of parametric phenomena , and our results are 
in good agreement with theoretical predictions. In 
particular , it may be used to study the influence of pump 
P waves on the instability’s threshold and amplitude of S 
waves. 
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I. INTRODUCTION 

The difference in speed of travel of P-waves and S-wave 
is vital to transmit energy of seismic wave. The P wave is 
a longitudinal wave or a compression wave. Force is 
applied in the direction that the wave is travelling. 
Ground or earth is pretty incompressible, so the energy is 
transferred pretty quickly. In S wave, the medium is 
displaced in a transverse (up and down - compared to the 
line of travel) way, and the medium must shear or "move 
away" from the material right next to it to cause the shear 
and transmit the wave [1,2]. 

On the other hand, parametric amplifiers and oscillators 
have been widely studied in electronics and optics. For 
instance, parametric amplification has been used to 
achieve low-noise amplification in electronic systems. 
This parametric instability is called the Faraday 
instability. Recently, Pritchett and Kim proposed a simple 
system to observe Faraday instabilities [3]. We propose 
here an alternative way to study parametric instabilities 
by doing an analog experiment that models the Mathieu 
equation [4]. 

A study can be made on the torsional-lateral motions of 
non-linear symmetrical structures subjected to lateral 
ground motion. The torsional and lateral response of a 
single mass symmetrical system subjected to sinusoidal 
ground motion can be investigated where non-linear 
coupling exists between the lateral and rotational 
motionslike S waves. For sinusoidal lateral response, the 


torsional motion equation can be cast in the form of a 
Mathieu equation. The likelihood of induced torsional 
response can be studied in terms of unstable regions in the 
parametric amplitude-frequency parameter space. The 
implication of this type of non-linear torsional-lateral 
coupling to the responses of real symmetrical structures 
subjected to actual earthquake ground motion can be 
simulated with an electronic model. 

Our proposed simulation is easy to understand 
conceptually and has several advantages, including fast 
data acquisition. It also allows students to explore the 
behavior of a driven oscillator and to understand the 
concepts of supercritical and subcritical bifurcations in 
earthquake phenomena. 

The experiments involve two control parameters: the 
forcing amplitude P and the forcing pulsation 00 e . 

Here, for wave P we make <D e = 2rad / S Varying these 

two parameters allows us to study the threshold of the 
instability for different driving frequencies and to explore 
the Mathieu extension of the bifurcation. We also study 
the nonlinear dependence of the oscillation amplitude on 

P. 

II. PARAMETRIC RESONANCE 

A system is subjected to a parametric forcing if one of its 
parameters is temporally modulated producing 
parametric instability which occurs when a tank 
containing a liquid is vertically vibrated: one then 
observes standing waves on the free surface. [5]. This 
parametric instability is called the Faraday instability. The 
physical situation is more complex because of the great 
number of degrees of freedom in the system leading to the 
generation of complex patterns on the surface. The study 
of parametric surface waves has led to a large number of 
theoretical and experimental studies [6]. In the simple 
case of an incompressible, irrotational, and inviscid fluid, 
Benjamin and Ursell [7], showed that, in the linear 

approximation, each mode £, k (of wave vector k) of the 
surface deformation is governed by a Mathieu equation. 
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+ G> 2 (l + Fcos(co e t)^ k =0 (1) 

at at 

which can be extended to rotational motion 

^ + :k^ + G3 2 (l + Pcos(co e t)L 1 =0 (2) 

dt 2 dt 0 e 0 

Where £, 0 is transversal This system leads to a canonical 

example of a parametric instability [5], 00 e is the external 

forcing pulsation, and P is directly related to the 
amplitude of the vibration acceleration relative to the 
acceleration of gravity produced by P- waves. The 
presence of a small viscous dissipation can be taken into 
account by including a phenomenological damping term 

X. 

In the undamped case X = 0 , when P — » 0 , the 
parametric resonance occurs when C0 e / <D 0 = 2 / n , 
where n is an integer. The most unstable oscillation 
corresponds to n = 1 , that is, C0 e / CO 0 = 2 , [8]. In the 
following, we are interested only in this last case. 

With X = 0 and CD 0 P = 8 Mathieu equation (2) is 

d 2 t 

— + (cOq + £ cos(2t))^ e = 0 , cOq > 0 (3) 

dt 

which is reversible, A typical question is: for which 
values of C0 Q and Sin ( CO 0 , 8 )-parameter space is the 

trivial solution = 0, d^ 0 / dt = 0 stable? 

We find that periodic solutions exist for n = 1 if: 
CD 0 =1 ±s/2 + O(s 2 ). In the case n - 2, periodic 

solutions exist if: OO 0 = 4 — S 2 / 48 + 0(s 4 ) , and 
co 2 =4+e 2 5/48 + 0(£ 4 ). 

The corresponding instability domains are called Floquet 
tongues, instability tongues or resonance tongues, see fig. 
1. 

On considering higher values of n , we have to calculate to 
a higher order of 8 . At n = 1 

the boundary curves are intersecting at positive angles at 
8 = 0 ,atft = 2(CO 0 =4) they are tangent; the order of 

tangency increases as n - 1 (contact of order n ), making 
instability domains more and more narrow with 
increasing resonance number n. 
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Fig.l : Floquet tongues of the Mathieu eq. (3); the 
instability domains are marked with /. 

III. THE MATHIEU EQUATION WITH VISCOUS 

DAMPING 

In real seismic applications there is always the presence 
of damping. We shall consider the effect of its simplest 
form, small viscous damping. Eq. (3) is extended by 
adding a linear damping term X : 

(4) 

We assume that the damping coefficient is small, 
X = 8K 0 , and we put CD 0 = n 2 — sP to 

apply the Poincare-Lindstedt method [9]. 

We find periodic solutions in the case n = 1 if: 

CO 2 = l ±^£ 2 / 4-^ 2 ( 5 ) 

Relation ( 5 ) corresponds with the curve of periodic 
solutions, which in ( CO 0 , 8 )-parameter 

space separates stable and unstable solutions. We observe 
the following phenomena. 

ifO<A,<e/2 , we have an instability domain which by 
damping has been lifted from the 

CO 0 -axis; also the width has shrunk. If -A,>s/2 the 

instability domain has vanished. For an 
illustration see fig. 2. 
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Fig. 2: Reduced schematic instability with the damping 
presence. 

Repeating the calculations for n > 2, we find no 
instability domains at all; damping of 0( 8 ) stabilizes the 
system for 8 small. To find an instability domain we 
have to decrease the damping, for instance if n = 2 we 

have to take • X = 8 2 K 0 (Figure 3). 



2 

Fig. 3: ( CD 0 , 8 ) space for different values of damping 


IV. AN OSCILLATORY ELECTRIC CIRCUIT 

Previous simulation can be verified experimentally 
through the following approach which can be used for a 
depth investigation of seismic S waves which can be 
attenuated or amplified. This proposed circuit and its 
dynamics can be approximately modeled with the 
Mathieu equation (4). 



Fig.4: Electronic analog of the Mathieu equation. 

We will propose an analog model of the parametric 
instability with an scaled frequency such that 

C0 e / CO 0 = 2 . In fact, the Mathieu equation, Eq. (4) 

above, leads to the amplitude in Eq. (5), which is 
analogous to the one derived in the Faraday problem in 
the limit of small viscosity. Our electronic device is 
shown in Fig 4. M is an AD633 multiplier; its output 

voltage, U 0 (t) is proportional to the product of the two 

input voltages, U 1 (t)xu 2 (t). C is a capacitor with a 

voltage-dependent capacity made with variable 
capacitance diodes (varicaps) BB909A. As mentioned in 
Ref. [10], these diodes introduce some nonlinearities in 
contrast to the circuit proposed in Ref. [11,12] where no 
nonlinear component is present. To avoid electromagnetic 
perturbations, the electronic oscillator can be enclosed in 
a metallic container which serves as a Faraday cage. 

The electronic system is forced by a periodic voltage 
U e (t) = U 0 COS(C0 e t) , C0 e / CO 0 = 2 . We can use 

Kirchhoff s voltage law and the property of the multiplier 
to find that the charge q of the capacitance C is governed 
by (see Fig. 4). The circuit LC have a time varying 
resistance R 

L^-^ + R — + — (l + ku n cos(2t)) = 0 (6) 

dt 2 dt C 0 

where k is the gain of the multiplier (k=l/10V) in our 
case. The capacitance C depends on the charge q, but the 
use of a pair of oppositely polarized varicaps leads to the 
symmetry q— »— qof the function C(q). Hence, we can 
write 

C = C 0 (7) 

+ K -7T + ®o(l + p e cos(2t))q = 0 (8) 

dt dt 

with X Q = R / L , cOq = 1 / LC 0 , 03qP 0 = oc>Qku 0 = 8 e . 
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““T + “T" + (®o + s e cos (2t:))q = 0 (9) 

dt dt 

Equation (9) is identical to expression for seismic waves 
and we can expect experimental results for the charge 
amplitude similar to those described for the P and S 
waves. 

We propose a 1-H inductance, and measured the 
resonance frequency of the RLC circuit to be 

oo 0 = 2n x 6.8 12kHz so that C G « 546pF . For 
CD e = 2 cD o , the threshold is simply determined by the 
relation P = X = R / (2L(D 0 ) . Thus, the global 
resistance of the circuit is R = 555Q, and the quality 
factor of the RLC circuit is QP = Lco 0 / R □ 77 . As 
mentioned below, we are interested in the evolution of the 
voltage U « q / C 0 , the analog of the amplitude £, 0 . We 

thus can utilize a digital synthesizer like HP8904A or 
equivalent and a signal analyzer HP35670A or equivalent 
to calculate, via a real-time averaged fast Fourier 

transform, the voltage amplitude | A| , defined by 

u(t) = Ae i( ° et/2 +cc (10) 

After interfacing the experimental setup using LAB VIEW 
or equivalent, we can directly record the control 

parameters U 0 with CD e = 2cO o , and the voltage 

amplitude | A| . Thus we can study the interaction of P-S 

waves [13], especially when the P-S interaction occurs in 
the metamaterial region. 

V. CONCLUSION 

We have proposed an experimental study of parametric 
resonance between P-waves and S-waves, which can be 
used to describe various nonlinear phenomena 
qualitatively and to obtain bifurcation diagrams 
quantitatively. It is based on an electronic circuit and is 
easy to design. We have shown that it is a good 
simulation of parametric phenomena, and our results can 
be in good agreement with theoretical predictions. In 
particular, it may be used to study the influence of pump 
P waves on the instability’s threshold and amplitude of S 
waves. 
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